Magnetic ground state of coupled edge-sharing CUO2 spin-chains 
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By means of density functional theory, we investigate the magnetic ground state of edge-sharing 
Cu02 spin-chains, as found in the (La,Ca,Sr)i4Cu24 04i system, for instance. Our data rely on spin- 
polarized electronic structure calculations including onsite interaction (LDA-I-U) and an effective 
model for the interchain coupling. Strong doping dependence of the magnetic order is characteristic 
for edge-sharing Cu02 spin-chains. We determine the ground state magnetic structure as function 
of the spin-chain filling and quantify the competing exchange interactions. 
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Doping of magnetic systems nowadays attracts great 
interest due to a large number of exotic effects. 3D anti- 
ferromagnetic order, for example, is easily destroyed by 
either quantum or thermodynamical fluctuations, and by 
any kind of impurity. The study of doped 2D antiferro- 
magnets fT] has been triggered by the phase diagrams of 
high-temperature superconductors. Upon doping, they 
show a transition from an antiferromagnetic to a super- 
conducting state. On the other hand, ID copper oxide 
compounds are an active field of research, too. From the 
theoretical point of view, the magnetic properties of un- 
doped two-leg CU2O3 ladders are explained in terms of a 
Heisenberg model with cyclic ring exchange 0. Doping 
effects are investigated in Refs. Antiferromagnetic 
CuO chains with 180° Cu-O-Cu bond angles, as found in 
Sr2Cu03 @, for example, can be understood by confor- 
mal field theory When a few Cu ions are replaced 
by non-magnetic ions, like Zn, the system is well de- 
scribed by the Heisenberg model with removed spins Q ■ 
Edge-sharing Cu02 chains with 90° Cu-O-Cu bond an- 
gles, on the contrary, are frustrated due to competition 
between ferromagnetic and antiferromagnetic exchange 
d, H, In this context, both charge and spin order of 
highly doped systems are not understood by far. 

The (Sr,Ca,La)i4Cu2404i system of isostructural ma- 
terials is prototypical for doping dependent charge and 
magnetic ordering in a partially filled spin-chain, i.e. a 
chain with less than one spin per site. The compounds 
have been subject of intensive research in recent years, 
mainly due to a rich phase diagram and a close relation 
to cuprate superconductors. The incommensurate crys- 
tal structures consist of planes of quasi ID Cu02 chains 
running along the crystallographic c-axis, stacked alter- 
nately with planes of two-leg CU2O3 ladders. The lattice 
constants of these subsystems satisfy lOcchain ~ 7ciaddcr 
[ll|, Materials with different Sr/Ca composition are 
isostructural with only minor modifications of the bond 
lengths and angles [ij, [ij, [l^ . A relative shift between 
adjacent chains parallel to the chain axis distinguishes Sr 
and Ca rich compounds. In Cai3.6Sro.4Cu2404i the shift 
amounts to half the intrachain Cu-Cu distance, while it 
is only 30% of this distance for Sri4Cu2404i. For each 
Sr/Ca composition the Cu valence is 2.25, where Cu^"*" 



ions {S — 1/2) accumulate on the ladders. In contrast, 
most of the 6 Cu'^"'" ions per unit cell, which form 5* = 
Zhang- Rice singlets (so-called holes), are located on the 
chains. Optical conductivity and x-ray absorption data 
suggest that substitution of Ca for Sr induces a transfer 
of holes back into the ladders [l^ [l3] ■ Consequently, an 
average Cu valence of 2.50 is assumed for the chain sites 
in Ca rich samples, giving rise to a half filled chain with 
one spin per two sites. On La substitution the intrinsic 
doping decreases gradually, reaching the undoped state 
with a nominal Cu valence of 2.00 for La6Ca8Cu2404i. 

By chain geometry, competing ferromagnetic nearest- 
neighbour and antiferromagnetic next-nearest neighbour 
coupling is expected for edge-sharing CUO2 spin-chains. 
The magnetic phase diagram comprises at least 3 phases 
upon doping, i.e. variation of the number of holes in the 
spin-chain. In the case of Sri4Cu2404i the Cu02 chains 
are non-magnetic with a spin gap of about 130 K, which 
can be explained in a dimer picture [l^ . Using a cluster 
calculation, Gellee and Lepetit [To'l have shown that the 
ground state is very complex, with trimer contributions. 
A dimer ground state is found by exact diagonalization of 
a Heisenberg model with nearest and next-nearest neigh- 
bour interaction ^2Q] . in accordance with band structure 
data [2l|. On the contrary, antiferromagnetic intrachain 
ordering has been reported for Ca rich samples [l^ [13] ■ 
Despite a huge number of investigations, most details of 
the antiferromagnetic order are still a matter of contro- 
versy. Even worse, in Sri3LaCu2404i spin dimers evolve 
instead of the antiferromagnetism expected for half filled 
chains ^231. Reflecting Cu-O-Cu bond angles of w 90°, 
ferromagnetism finally is realized in slightly doped spin- 
chains of La rich systems [l^, , which are captured by 
a 2D classical Heisenberg model with defects [26| . 

Spin-chain systems are strongly affected by both sub- 
tleties of the crystal structure and the electron-electron 
interaction. In the case of half filled spin-chains [27| the 
ground state magnetic structure has been established in 
[23] ■ However, the dependence of the exchange coupling 
on the filling n (where n = number of spins / number of 
sites) is an open question, since the coupling over a hole 
competes with the coupling over an interstitial spin. In 
order to obtain a comprehensive picture of doped spin- 
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TABLE I; Total energy gain (per Cu atom) with respect to a 
non-magnetic calculation, for both staggered and symmetri- 
cally aligned spin-chains. The fm, afm, 4kF/afm, and cluster 
spin patterns are compared to each other [2^]. The filling of 
the spin-chain is n = 1/2. 



chains, we thus clarify the details of the transition from 
antiferromagnetic to ferromagnetic coupling on reducing 
the hole count. To be more specific, we study the onset 
of ferromagnetism at n = 2/3 and the full spin chain at 
n — 1 for typical Cu-Cu bond lengths of 2.75 A (nearest 
neighbour) and 5.8 A (next-nearest neighbour). 

The following results are based on density functional 
theory and the generalized gradient approximation. We 
apply the WIEN2k program package, a state-of-the-art 
full-potential code with a mixed lapw and apw-|-lo basis 
[29( 1 . For all the spin configurations under investigation, 
the charge density is represented by at least 70,000 plane 
waves and the Brillouin zone integration uses a k-mesh 
comprising at least 18 k-points in the irreducible wedge. 
While Cu 3p and O 2s orbitals are treated as semi-core 
states, the valence states contain Cu 3c?, 4s, 4p and O 
2p, 3s orbitals. To achieve a mixed valence in a unit cell 
with both Cu^+ and Cu^+ ions, we start the calculation 
from the following setup: One electron is added to each 
oxygen, thus preventing charge transfer off the copper, 
and restricted to a Watson sphere of radius 1.3 A. More- 
over, a 3(i^° configuration is used for the Cu'^+ ions and 
a id^'^As^ configuration for the Cu^+ ions. 

From a structural point of view, the spin-chain mate- 
rials (La,Ca,Sr)i4Cu2404i consist of 3 subsystems: The 




ff/f 



ff/a 



fa/f 



fa/a 



4) 4, 4) 
4-4,4' 
4-4,4' 
4-4,4- 
4-4,4- 
4-4,4) 

3f3a/f^C()^ 
3f3a/a^ (() ^ 
2fla/a^ (() ^ 



f/a 



4)°o*(t)4'4)°o 

4.°o*4)*4)°o 
4)°o*4)*4)°o 
4)°o*4)*4)°o 



4-*4- 



FIG. 1: Density of states (DOS) for the chain subsystem of 
Sri4Cu2404i. [33] 



FIG. 2: Nomenclature of spin patterns on coupled chains. 



Cu02 spin-chains, CU2O3 ladders, and interstitial elec- 
tron donor ions. As interaction between the subsystems 
is very small, an effective model of coupled spin-chains, 
capturing structural prerequisites of the exchange inter- 
action, may be used to investigate the compounds [2l| . 
The density of states (DOS) of the Cu02 chains shows a 
structure of width Ki 5.2 eV at binding energies larger 
than 2 eV and a second structure of width 1 eV right 
at the Fermi level. As expected for transition metal 
oxides [13, HH [111, the states are due to bonding be- 
tween Cu 3dxz and O 2px/pz orbitals. Characteristic Cu 
idxz bands at the Fermi energy dominate the electronic 
and magnetic properties [1^. A tight-binding fit reveals 
strong next-nearest neighbour hopping, which is rather 
independent of the interchain coupling. Nearest neigh- 
bour hopping, in contrast, can be suppressed by relative 
shifts between the Cu02 chains. In that case the trans- 
verse hopping, a measure of the 2D coupling, increases 

For n — 1/2, the local density approximation (LDA), 
not accounting for the onsite interaction, yields a ferro- 
magnetic ground state for staggered spin-chains [28^. In 
contrast, antiferromagnetic spin order and 4kF-periodic 
charge order, complemented by antiferromagnetic inter- 
chain coupling, are observed for symmetric spin-chains, 
compare the data in Table [T| In this case, the ferromag- 
netic solution is even unstable. Including local electron- 
electron interaction at the correlated Cu id orbitals, an 
LDA-I-U treatment reveals an antiferromagnetic ground 
state in both arrangements, see Table [T| A relative shift 
between the Cu02 chains promotes antiferromagnetism 
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TABLE 11: Total energy gain (per Cu atom) with respect to a noii- magnetic calculation for the spin-chain filling factor n — 2/3. 
Several spin patterns are compared to each other, see fig. [2] for the nomenclature. 
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TABLE IIL Total energy gain (per Cu atom) with respect to 
a non-magnetic calculation for a full spin-chain (n = 1). The 
f/f, f/a, and a spin patterns are compared to each other. 

more effectively than a local interaction. In total, next- 
nearest neighbour coupUng over a hole site controls the 
magnetism in a half filled spin-chain. The combination 
of electron-electron repulsion and structural frustration 
therefore gives rise to a 4kp spin pattern. Due to a weak 
interchain coupling of about 1 mRyd, the magnetic next- 
nearest neighbour interaction energy is 54 ± 2 mRyd per 
Cu atom within the chain. 

Turning to a lower hole doping, we start with the fill- 
ing n = 2/3 and afterwards consider the full spin-chain, 
given hy n = 1. For n — 2/3, we have 4 Cu^+ ions per 6 
lattice sites on the coupled chains. The supercell enter- 
ing the band structure calculation therefore contains 12 
inequivalent Cu sites. We set up the start configuration 
as follows: Within a chain, a ferro- or antiferromagnetic 
two spin cluster is coupled to the next such cluster over 
a hole site. Experimental data rather point at antiferro- 
magnetic than ferromagnetic coupling between adjacent 
chains [1^. However, as the details of the exchange in- 
teraction are not known, we study both possibilities for 
comparison. This leads to altogether five configurations, 
which we systematically call ff/f, fa/f, ff/a, fa/a and a. 
For clarity, these configurations are depicted in fig. [21 In 
each case, we calculate the total energy and compare it 
to the findings of a spin-degenerate calculation, yielding 
the energy gain due to the exchange interaction. 

The fully ferro- (ff/f) and antiferromagnetic (a) cluster 
spin configurations in the following serve as a reference 
frame for the discussion. Without local electron-electron 
interaction, only these two patterns are found to be sta- 
ble. The magnetic moments are « 0.35/iB on the Cu^+- 
sites, where moments connected to two holes are slightly 
larger than those connected to a single hole. Ferromag- 
netic clusters with antiferromagnetic intrachain coupling 
(fa/f and fa/a) transform into modified clusters without 
intermediate holes, therefore into a more homogeneous 
charge distribution. According to the initial interaction, 
we obtain a 3f3a/f or 3f3a/a spin pattern, as depicted in 
fig. O Ferromagnetic intrachain coupling in combination 



with antiferromagnetic interchain coupling (ff/a) results 
in a homogeneous spin-chain (fa). For the spin patterns 
under consideration, table [IT] compares the total energy 
gain due to the magnetic exchange interaction. 

When we account for a local electron-electron interac- 
tion U = 0.3 Ryd by means of the LDA-I-U method, the 
spin patterns ff/f, 3f3a/f, 3f3a/a and f/a remain stable, 
where the magnetic moments increase slightly to about 
0.5 /Xfi. Only the fully antiferromagnetic cluster configu- 
ration develops a more homogeneous charge distribution 
given by the spin pattern 2fla/a, as shown in fig. [51 We 
obtain the largest total energy gain (with respect to the 
non-magnetic state) for the f/a spin pattern, amounting 
to 32 mRyd per Cu site. The magnetic ground state of 
edge-sharing Cu02 chains with filling n = 2/3 therefore 
is given by fully polarized ferromagnetic spin-chains and 
antiferromagnetic interchain coupling. Remarkably, the 
electronic correlations do not alter the magnetic ground 
state, see the U = ORyd data in table [Til As compared 
to the coupling in a half filled spin-chain, ferromagnetic 
nearest neighbour interaction here becomes dominant as 
the electron-electron repulsion supports a homogeneous 
charge distribution. Whereas we find the onset of ferro- 
magnetism already for n — 2/3, experiments on La-rich 
systems suggest a nominal filling of about n = 0.7. We 
ascribe this discrepancy to insufficient knowledge about 
the effective chain hole count in the real materials. By 
the data given in table [Til we conclude that the energy 
gain due to the interchain magnetic coupling is at most 
1 mRyd per lattice site. In addition, comparison of the 
f/a to the 2fla/a data allows us to estimate an energy 
difference of 1.3 mRyd between the ferro- and antiferro- 
magnetic coupling of neighbouring spins. The combina- 
tion of the ff/ f data with the half filled spin chain results 
in an intrachain nearest neighbour exchange interaction 
energy of 27 ±2 mRyd per Cu atom. Therefore, the next- 
nearest neighbour coupling is found to be almost exactly 
twice as strong as the nearest neighbour coupling. 

For the full spin-chain (n = 1), we address ferro- and 
antiferromagnetically coupled Cu02 chains with a ferro- 
magnetic intrachain coupling (f/f and f/a), and a frus- 
trated antiferromagnet (a). For all three configurations 
the band structure calculation converges with magnetic 
moments of about 0.4 /ib- The antiferromagnet has the 
lowest energy, both for U = Ryd and a moderate local 
electron-electron interaction (of U = 0.3 Ryd), see table 
ITTTl In the latter case, the magnetic moments increase to 
0.6 fiB for the antiferromagnet. In contrast to the doped 
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compounds, a local interaction now counteracts the spin 
ordering. By an energy difference of less than 0.2 mRyd 
per lattice site, ferro- and antiferromagnetic interaction 
between nearest neighbour spins result in almost degen- 
erate electronic states. On the contrary, the energy gain 
due to the interchain coupling is larger than determined 
for the spin-chain fillings n = 1/2 and n = 2/3. 

In conclusion, we have investigated the magnetic cou- 
pling in edge-sharing Cu02 chains with sa 90° Cu-O-Cu 
bond angles, at doping levels of n = 1/2, n = 2/3, and 
n = 1. By means of LDA-I-U electronic structure calcu- 
lations, we have studied the total energy for various spin 
patterns. Our data shows that accounting for both the 
chain geometry and the electron-electron interaction [s^l 
is essential for obtaining a comprehensive picture of the 
magnetism, though electronic correlation does not alter 
the magnetic ground state. For symmetric Cu02 chains 
and filling n = 1/2 the latter is antiferromagnetic with a 
4kF charge order along the chain axis. Adjacent chains 
are coupled antiferromagnetically, too. This result par- 
ticularly applies to the coupled Cu02 chains of Ca-rich 
spin-chain materials (Sr,Ca)i4Cu2404i. La-rich systems 



with spin-chain filling n = 2/3 develope a ferromagnetic 
intrachain coupling since Coulomb repulsion supports a 
homogeneous distribution of the charge. The interchain 
coupling remains antiferromagnetic. When the filling is 
complete {n — 1), finally a frustrated antiferromagnet is 
found, which fact reflects strong electronic correlations. 
Under the assumption that modelling in terms of a Hub- 
bard model is appropriate, antiferromagnetism supports 
intrachain hopping, which is excluded for ferromagnetic 
spin-chains. By Nagaoka's theorem [s^, a single hole is 
expected to turn the system into a ferromagnetic state, 
in agreement with results for La5Ca9Cu2404i. However, 
the experimental observation of antiferromagnetic spin- 
chains in undoped La6Ca8Cu2404i is missing so far. 
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